Positive and negative Hanbury-Brown and Twiss correlations in normal 

metal— superconducting devices 
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In the light of the recent analogs of the Hanbury-Brown and Twiss experiments 111 in mesoscopic 
beam splitters, negative current noise correlations are recalled to be the consequence of an exclusion 
principle. Here, positive (bosonic) correlations are shown to exist in a fermionic system, composed 
of a superconductor connected to two normal reservoirs. In the Andreev regime, the correlations 
can either be positive or negative depending on the reflection coefficient of the beam splitter. For 
biases beyond the gap, the transmission of quasiparticles favors fermionic correlations. The presence 
of disorder enhances positive noise correlations. Potential experimental applications are discussed. 
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In condensed matter systems, correlations effects be- 
tween carriers exist either because particles interact with 
each other, or alternatively because the observable one 
considers involves a measurement on more than one par- 
ticle. The characterization of current fluctuations in time 
constitutes a central issue in quantum transport. Noise 
measurements have been used to detect the fractional 
charge of the excitations in the quantum Hall effect j|,|J . 
More recently, a fermion analog of the Hanbury-Brown 
and Twiss experiment J2[] was achieved [|l| with meso- 
scopic devices, obtaining a clear signature of the nega- 
tive correlations expected from the Pauli principle. Here 
we recall rapidly the ingredients which are necessary for 
negative correlations and we propose a Hanbury-Brown 
and Twiss experiment for a fermionic system were both 
negative and positive (bosonic) noise correlations can be 
detected. 

The system which is proposed (see inset of Fig. [j]) con- 
sists of a junction, or beam splitter, connected by electron 
channels to reservoirs, which is similar to that of Ref. |]6), 
except that the injecting reservoir is a superconductor. 
Because of the proximity effect at the interface between 
the superconductor and the normal region, electrons and 
holes behave like Cooper pairs provided there is enough 
mixing between them. While Cooper pairs are not bosons 
strictly speaking, an arbitrary number of theses can exist 
in the same momentum state, which opens the possibil- 
ity for bosonic correlations. Bosonic behavior in electron 
systems has been previously discussed, for instance for 
excitons in coupled quantum wells pi, where the pos- 
sibility of observing Bose condensation is in debate. On 
the other hand, it may be possible to detect negative cor- 
relations in adequately prepared photonic systems |S| . 

Negative noise correlations in branched electron cir- 
cuits are the consequence of an exclusion principle, which 
exists for fermions (Pauli principle) or even for particles 
which obey Haldane's exclusion statistics H. In these 
two situations, MMO, ignoring thermal effects, low fre- 
quency current noise in a two-probe device is suppressed 
by a factor (1 — T), with T the transmission probability. 



Consider first a circuit with three leads with correspond- 
ing currents and noises Ii and Si labeled i = 1,2,3, as 
depicted in inset of Fig. |l| (but ignoring region 4, which 
we take later to be a superconductor) . Assume that par- 
ticles (fermions or exclusion particles) with charge q are 
injected from 3 with a chemical potential /13 while 1 and 
2 are kept at the same chemical potential fi± . Noise cor- 
relations between 1 and 2 can be computed by invoking 
that the fluctuations in 3 M equal that of a composite 
lead (1 + 2): at u — 0, S3 — 5 , ( 1+2 )- The definition of the 
noise correlations between 1 and 2 is: 
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with Sli the fluctuation around the average current in i. 
The correlations S12 are obtained by subtracting the indi- 
vidual noise of 2 and 3 from 1: S 12 = [S3 - (Si + S 2 )}/2. 
For electrons and exclusion particles, a multi-terminal 
noise formula M gives: 
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were S(i+2),3 is the (2 x 1) transmission matrix between 
3 and (1 + 2), which is a submatrix of the scattering ma- 
trix with elements Sij describing the junction, 1 is the 
two-dimensional identity matrix, and g is the exclusion 
parameter (g = 1 for fermions). Using current conserva- 
tion, one obtains negative correlations for both fermion 
and exclusion particles: 
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Minimal negative correlations S12/S1 — —1 are obtained 
for a reflectionless, symmetric junction. 

Positive correlations in systems where the inject- 
ing lead is a superconductor are now addressed. 
The scattering approach to quantum transport in the 
presence of normal-superconductor (NS) boundaries 
is available |l2]-[il|] so the basic steps are reviewed 



briefly. The fermion operators which enter the cur- 
rent operator are given in terms of the quasipar- 
ticle states using the Bogolubov transformation Jig ] 

Vvfc) = En [uni^Cna ~ (TV* (x)c^ _ a j , were 4<, (Cntr) 

are quasiparticle creation (annihilation) operators, n = 
(i, a, -E) contains information on the reservoir (i) from 
which the particle (a — e,h) is incident with energy E 
and a labels the spin. The contraction of these two op- 
erators gives the distribution function of the particles in- 
jected from each reservoir, which for a potential bias V 
are: fi e = f(E — eV) for electrons incoming from i, simi- 
larly fih = f(E + eV) for holes, and /i ]Q = f(E) for both 
types of quasiparticles injected from the superconductor 
(/ is the Fermi-Dirac distribution). Here, eV > means 
that electrons are injected from regions 1 and 2. Invoking 
electron-hole symmetry, the (anti)correlations of holes 
are effectively studied. u n (x) and v n (x) are the solutions 
of the Bogolubov-de Gennes equations which contain 
the relevant information on the reflection/transmission 
of electrons and holes (and their quasiparticle analogs) at 
the NS interface. The current operator allows to derive 
a general expression for the zero frequency noise corre- 
lations between normal terminals i and j U4R& which 
constitutes our starting point: 
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where current matrix elements are defined by Ai a jp = 
Ujpd x u* a - u* a d x u lf3 and B iaj/3 = v* jl3 d x v ia - v ia d x v* . 
The electron and hole wave functions describing scatter- 
ing states a (particle) and i (lead) are expressed in terms 
of the elements Sij a p of the S-matrix which describes the 
whole NS ensemble: 

(5a) 
(5b) 

where Xj denotes the position in normal lead j and 
fc± ( v ±) are the usual momenta (velocities) of the two 
branches. Sij(0) has been shown to have no definite sign 
in four-terminal noise measurements E3]. 

Specializing now to the NS junction connected to a 
beam splitter (inset of Fig. Q), 6 x 6 matrix elements 
are sufficient to describe all scattering processes. At zero 
temperature, the noise correlations between the two nor- 
mal reservoirs simplify to: 
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where the subscript 4 denotes the superconducting lead. 
The first term represents normal and Andreev reflec- 
tion processes |T7J] , while the second term invokes the 
transmission of quasiparticles through the NS boundary. 
It was noted previously Jig ] that in the pure Andreev 
regime the noise correlations vanish when the junction 
contains no disorder: electron (holes) incoming from 1 
and 2 are simply converted into holes (electrons) after 
bouncing off the NS interface. The central issue, whether 
disorder can induce changes in the sign of the correla- 
tions, is now addressed. 
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FIG. 1. Noise correlation between the two normal reser- 
voirs of the device (inset), as a function of the transmission 
probability of the beam splitter, showing both positive and 
negative correlations. Inset: the device consists of a super- 
conductor (4)-normal (3) interface which is connected by a 
beam splitter (shaded triangle) to reservoirs (1) and (2) 



Consider first the pure Andreev regime, were eV <C A, 
the superconducting gap, for which a simple model for a 
disordered NS junction 1 13] is readily available. The junc- 
tion is composed of four distinct regions (see inset Fig. 
1). The interface between 3 (normal) and 4 (supercon- 
ductor) exhibits only Andreev reflection, with scattering 
amplitude for electrons into holes ta = 7exp(— i(f>) (the 
phase of 7 = exp[— i arccos(i?/A)] is the Andreev phase 
and (j) is t ne phase of the superconductor). Next, 3 is 
connected to two reservoirs 1 and 2 by a beam splitter 
which is parameterized by a single parameter < e < 1/2 
identical to that of Ref. |pj: the splitter is symmetric, its 
scattering matrix coefficients are real, and transmission 
between 3 and the reservoirs is maximal when e = 1/2, 
and vanishes at e = 0. Electrons and holes undergo mul- 
tiple reflections in central 3 and the scattering matrix 
coefficient Sij a p of the whole device are computed using 
the analogy with a Fabry-Perot interferometer. 
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were x = yl — 2e and the remaining coefficients of the 
S-matrix are found using time reversal symmetry. 

Next one proceeds with the standard approximation 
7 ~ — i which applies for low biases in order to perform 
the energy integrals in Eq. (O): 
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The noise correlations vanish at e = 0, when conductors 
1 and 2 constitute a two-terminal device decoupled from 
the superconductor, and in addition, S±2 vanishes when 
e = \/2 — 1. A plot of S12 as a function of the beam 
splitter transmission (Fig. [j]) indicates that indeed, the 
correlations are positive (bosonic) for < e < \/2 — 1 
and negative (fermionic) for v2— 1 < e < 1/2. At max- 
imal transmission into the normal reservoirs (e = 1/2), 
the correlations normalized to the noise in 1 (or 2) give 
the negative minimal value: electrons and holes do not 
interfere and propagate independently into the normal 
reservoirs. It is then expected to obtain the signature of 
a purely fermionic system. When the transmission e is de- 
creased, Cooper pairs can leak in region 3 [Q because of 
multiple Andreev processes. Further reducing the beam 
splitter transmission allows to balance the contribution 
of Cooper pairs with that of normal particles. Eq. (H) 
predicts maximal (positive) correlations at e = 1/3: a 
compromise between a high density of Cooper pairs and 
weak transmission. 

The model described above may not be convincing 
enough, as an ideal Andreev interface was assumed. 
Moreover, it does not allow to generalize the results to 
the case where quasiparticles in the superconductor con- 
tribute to the current. Quasiparticles have fermionic 
statistics, so their presence is expected to cancel the pos- 
itive contribution of Cooper pairs leaking on the normal 
side. In particular, in the limit were eV 3> A, one should 
recover fermionic correlations. 

These issues bear similarities with a recent discussion 
of singularities in the finite frequency noise of NS junc- 
tions |lj|: below the gap, a singularity exists at the 
Josephson frequency, while above gap there appear ad- 
ditional features at eV ± A associated with electron and 
hole-like quasiparticles, which give single-particle behav- 
ior in the limit eV 3> A. However, finite frequency noise 
probes the charge (e or 2e) of the carriers, while here one 
is probing the statistics of the effective carriers in the 
junction. 

The energy dependence of the scattering coefficients is 
therefore needed to describe the correlations away from 
the pure Andreev regime. The numerical calculations 
which follow are performed using the BTK model p0[, 



where the disordered interface between regions 3 and 4 
(inset of Fig. fl\) can be characterized by a small number 
of parameters: the pair amplitude is assumed to be a step 
function A(x) = A0(x) and a delta function potential 
barrier V(x) — HvpZS(x) is imposed, where vf the Fermi 
velocity and Z > 1 (Z < 1) for strong (weak) disorder. 
The beam splitter is taken to be similar to the previous 
calculation M , assuming that the reflection/transmission 
of electrons does not depend significantly on the incom- 
ing energy. 
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FIG. 2. Noise correlations using an NS boundary mod- 
eled by BTK for weak disorder, Z — 0.1: from top to bottom 
eV/A — 0.5,0.95, 1.2, 1.8. Inset: proposed device for the ob- 
servation of positive/negative correlations; at the boundary of 
a superconductor (S), two point contacts (Pi) and (P2) are 
connected to a semitransparent mirror (M) 

Consider the case of weak disorder, Z = 0.1 (Fig. 2). 
At weak biases, good agreement is found with the previ- 
ous analytical results displayed in Fig. W, except that for 
a fully transmitting splitter, the ratio of the correlations 
divided by the noise in region 1 does not reach the ex- 
tremal value —1: an early signature of disorder. When 
the bias is further increased but is kept below the gap, 
the phase accumulated in Andreev processes by electrons 
and holes with various energies is spread over the inter- 
val [0,7r/2]: positive correlations are weaker, but they 
survive at low beam splitter transmission. Further in- 
creasing the voltage beyond the gap destroys completely 
the bosonic signature of the noise. 

A strikingly different behavior is obtained for inter- 
mediate disorder at Z = 1 (Fig. M). First, for weak 
biases, the noise correlations remain positive over the 
whole range of e, with a maximum located at e ~ 0.43, 
which is close to the case of a reflectionless splitter. This 
maximum becomes a local minimum for higher biases, 
where positive correlations remain quite robust never- 
theless. Just below the gap (eV = 0.95), correlations 
oscillate between the positive and the negative sign, but 
further increasing the bias eventually favors a fermionic 



behavior. Calculations for larger values of Z confirm the 
tendency of the system towards dominant positive cor- 
relations at low biases with S , i 2 (e)/S'i(l/2) > 1 over a 
wide range of e (not shown) . The phenomenon of positive 
correlations in fcrmionic systems with a superconducting 
injector is thus enhanced by disorder at the NS boundary. 
Nevertheless, for strong disorder, the absolute magnitude 
of S\ and S\i becomes rather small, which limits the pos- 
sibility of an experimental check in this regime. 



and holes are correlated. This can be viewed as a fi- 
nite density of Cooper pairs which behave like bosons. 
The presence of disorder allows in some cases to enhance 
the appearance of bosonic correlations. Similar studies 
could be envisioned in the Fractional Quantum Hall Ef- 
fect (FQHE) where the collective excitations of the cor- 
related electron fluid have unconventional statistics M . 

Discussions with the late R. Landauer, with D.C. Glat- 
tli and with M. Devoret are greatfully acknowledged. 
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FIG. 3. Noise correlations using an NS boundary modeled 
by BTK for intermediate disorder: Z — 1 (same biases as in 
Fig. 2) 

An interesting feature of the present results is the fact 
that both positive and negative correlations are achieved 
in the same system. A suggestion for this device is de- 
picted in the inset of Fig. 2. Assume that a high mobility 
two dimensional electron gas has a rather clean inter- 
face with a superconductor pi| . A first point contact 
close to the interface selects a maximally occupied elec- 
tron channel, which is incident on a semi-transparent 
mirror similar to the one used in the Hanbury-Brown 
and Twiss fermion analogs H. A second point contact 
located in front of the mirror, allows to modulate the re- 
flection of the splitter in order to monitor both bosonic 
and fermionic noise correlations. In addition, by choosing 
a superconductor with a relatively small gap, one could 
observe the dependence of the correlations on the voltage 
bias without encountering heating effects in the normal 
metal. 

Hanbury-Brown and Twiss type experiments may be- 
come a useful tool to study statistical effects in meso- 
scopic devices. Here, noise correlations have been shown 
to have either a positive or a negative sign in the same 
system. Close to the boundary, a fraction of electrons 



[1 
[2. 
[3 
[4 
[5 
[6 

[r 

[8 

[9 

[10 

[11 
[12. 

[13 

[14 
[15 

[16 

[17 

[is; 

[19 

[20 
[21 



M. Henny et al., Science 284, 296 (1999); W. Oliver et 

al., ibid, 299 (1999). 

R. Hanbury-Brown and Q. R. Twiss, Nature 177, 27 

(1956). 

L. Saminadayar, D. C. Glattli, Y. Jin, and B. Etienne, 

Phys. Rev. Lett. 79, 2526 (1997). 

R. de-Picciotto, M. Reznikov, M. Heiblum, V. Umansky, 

G Bunin, and D. Mahalu, Nature 389, 162 (1997). 

Y. Gefen, Y. Imry and R. Landauer, Phys. Rev. Lett. 52, 

139 (1984). 

T. Martin and R. Landauer, Phys. Rev. B 45, 1742 

(1992). 

V. Negoita, D. W. Snoke and K. Eberl, Cond. Mat. 

9901088 (and references therein). 

R. Liu, private communication. 

F.D.M. Haldane, Phys. Rev. Lett. 67, 937 (1991). 

GB. Lesovik, JETP Lett. 49, 592 (1989); M. Biittiker, 

Phys. Rev. Lett. 65, 2901 (1990). 

S. Isakov, T. Martin and S. Ouvry, Cond. Mat. 9811391. 

C.W.J. Beenakker, in Mesoscopic Quantum Physics, eds. 

E. Akkermans et al., p. 279 (Les Houches LXI, North 

Holland 1995). 

B. A. Muzykantskii and D. E. Khmelnitskii, Phys. Rev. 

B 50, 3982 (1994). 

M. P. Anantram and S. Datta, Phys. Rev. B 53, 16 390 

(1996). 

P.G de Gennes, Superconductivity of Metals and Alloys, 

(Addison Wesley, 1966, 1989). 

T. Martin, Phys. Lett. A 220, 137 (1996). 

A.F. Andreev, J. Exp. Theor. Phys. 46, 1823 (1964) [Sov. 

Phys. JETP 19, 1228 (1964)]. 

See for example, A. A. Abrikosov, Fundamentals of the 

Theory of Metals, (North-Holland, 1988) 

G Lesovik, T. Martin and J. Torres, Cond. Mat. 9902278; 

J. Torres, G Lesovik and T. Martin, [in preparation, 

1999). 

GE. Blonder, M. Tinkham and T.M. Klapwijk, Phys. 

Rev. B 25, 4515 (1982). 

A. Dimoulas et al, Phys. Rev. Lett. 74, 602 (1995). 



